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1 .  Introduction 

The  Global  Grid  communications  environment  is  designed  to  provide  a  high-data-rate  network  that 
supports  the  full  complement  of  DoD  communications  needs.  Satellite  links  are  an  integral  part  of 
the  Global  Grid  concept  which  extends  the  connectivity  to  individual  mobile  units  and  jointly 
operating  forces  at  isolated  locations.  However,  satellite  channels  are  vulnerable  to  jamming,  in 
particular  when  a  high  data  rate  is  required.  To  overcome  the  jamming  problem,  the  traditional 
approach  is  to  design  modulation  techniques  such  that,  on  the  average,  the  spectrum  of  the  signal 
looks  like  the  spectrum  of  a  white  noise  process.  This  together  with  channel  coding  provides 
immunity  against  jamming.  The  main  problem  with  this  approach  is  that  it  requires  a  much  larger 
spectrum  than  that  of  the  data  sequence;  hence,  the  available  spectrum  on  the  link  is  used  mostly 
for  jamming  protection  and  not  for  data  communications. 

The  traditional  modulation  techniques  employed  on  satellite  communication  links  are  based  on 
various  properties  of  the  second-order  spectrum  of  the  modulation  wave  form.  In  this  study  we 
have  explored  the  utility  and  performance  of  a  new  family  of  modulation  schemes  that  exploit  the 
properties  of  the  higher-order  cumulant  sequences  and  associated  poly  spectra  of  the  waveform.  In 
particular,  we  have  investigated  an  approach  in  which  the  third-order  polyspectrum,  which  is 
generally  referred  to  as  the  bispectrum,  is  modulated.  To  determine  the  performance  characteristics 
of  this  bispectral  modulation  scheme,  we  have  considered  two  different  detector  structures,  for 
which  we  have  completed  both  theoretical  and  simulated  performance  analyses.  We  have  studied 
the  performance  of  each  detector  structure  in  both  an  additive  white  Gaussian  noise  environment 
and  a  partial-band  jamming  environment,  and  we  have  demonstrated  that  bispectral  modulation 
techniques  can  significantly  outperform  traditional  anti-jamming  modulation  schemes,  such  as 
frequency-hopped,  binary  frequency-shift  keying,  in  the  presence  of  partial-band  jamming.  Our 
results  also  indicate  that,  even  on  an  additive  white  Gaussian  noise  channel,  bispectral  modulation 
can  be  used  in  parallel  with  conventional  noncoherent  modulation  techniques  to  increase  the  data 
rate  on  the  channel  or  to  reduce  the  sensitivity  to  synchronization  errors  between  the  transmitter  and 
receiver. 

This  report  is  organized  as  follow.  In  Section  2,  we  present  a  brief  introduction  to  polyspectra  and 
one  particular  polyspectrum  estimation  technique.  In  Section  3,  we  discuss  in  detail  the  proposed 
modulation  scheme  and  a  simple  technique  for  implementing  it.  The  results  of  our  performance 
analyses  are  presented  in  Sections  4  and  5,  and  in  Section  6  we  present  some  conclusions  about  the 
research  and  suggestions  for  further  study.  Considerably  more  detail  regarding  the  advantages  and 
direction  of  future  research  will  be  presented  in  our  Phase  II  proposal,  which  will  be  submitted 
within  30  days. 


2 .  Background  and  Definitions 

In  this  section,  we  give  some  necessary  background  information  on  the  properties  of  higher-order 
moments  and  polyspectra  and  a  brief  discussion  of  the  polyspectrum  estimation  technique 
employed  in  this  study. 


Let  {x(n)}  be  a  strictly  stationary  real  random  process.  We  assume,  without  loss  of  generality, 
that  (jc(n)}  has  mean  zero.  The  &th- order  cumulant  sequence  for  such  a  process  is  defined  in  terms 
of  the  first  k  moments  of  the  process.  For  example,  the  second-,  third-,  and  fourth-order  cumulant 
sequences  are  given  by: 
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C2(n  i )  =  E{x(&).x(&  +  nx  )}, 

C?>{n\’nl)  ~  E ^x(k)x[k  +  rii)x(k  +  n2)}, 

C4  (rtj  ,n2,n2)  =  E{x(&).x(&  +  nx  )x{k  +  n2  )x[k  +  n2 )}  -  E{x(A:).x(fc  +  n\  )}e{x(&  +  n2  )x[k  +  n3 )} 

-E|jc(/:)a:(A:  +  n2  )}E  {*(/:  +  ri\  )x(k  +  n3 )}  -  E  j  x(k  )*(&  +  «3  )}E  jx:(&  +  n,\  )x(k  +  n2 )}, 

where  the  symbol  “E”  represents  the  expectation  operator,  and  k,  n\,  n2,  and  n3  are  arbitrary 
integers.  Clearly  {C2(«i)}  is  just  the  covariance  sequence  of  the  process  and  {Ci(n\,n2)}  is  the 
third  moment  sequence. 

The  Mi- order  polyspectrum  of  the  process  5^(0))  is  the  multidimensional  Fourier  transform  of  the 
corresponding  cumulant  sequence;  that  is, 

^(co)  =  ^(©1,---,toJt_1)=  f)  •••  2cfc(n1»*“'n*-1)g"^Clh"1+','+<ft‘"lB‘"^ 

«i=- 00  «t-i=-°° 

For  k  =  2,  the  polyspectrum  is  clearly  just  the  power  spectrum  of  the  process.  For  k  =  3  and  k  =  4 
(the  other  two  cases  of  most  practical  interest),  the  polyspectra  are  usually  referred  to  as  the 
bispectrum  and  trispectrum ,  respectively.  Due  to  the  symmetries  arid  periodicities  of  the 

polyspectrum,  it  is  completely  defined  by  its  values  on  a  nonunique  principal  domain  Dk  c  Rk  l. 
For  example,  the  principal  domain  of  the  bispectrum  is  generally  taken  to  be: 

Z>5  =  {(g)!,co2):0  <  tOj  ^  k,  0  <  oo2  -  minCo)!,  271-2(0!)}. 

This  region  is  depicted  in  Figure  1 . 
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Figure  1.  Principal  Domain  of  the  Bispectrum 

Polyspectra  have  many  interesting  properties,  but  the  following  three  are  perhaps  the  most  useful: 

1 .  If  {*(«)}  and  {y(n)}  are  independent  processes  with  polyspectra  S  xk  and  S  {,  then 

= st + si 

2 .  If  {jc in) }  is  a  Gaussian  process,  then  Sxk  =  0  for  k  >  2. 

3 .  If  {y(n) }  is  produced  as  the  output  of  a  linear  system  with  frequency  response  H((d)  driven 
by  a  process  {*(«)},  then 

Si (cp)  =  H{ co,  )H{ co2  )•  •  •  H{ak_x  )H{-v>x  -  to2 - (%_!  )Sl  (cp) . 

In  this  report,  we  will  be  concerned  exclusively  with  the  bispectrum,  which  we  will  denote  simply 
as  6(0)1,02).  The  following  examples  provide  some  simple  illustrations  of  the  bispectrum  and  its 
properties. 

Example  1  -  Linear  processes.  If  {£(«)}  is  an  i.i.d.  sequence  of  random  variables  such  that 

M-l 

then  Be(co1,to2)  -  1°  general,  if  {x(fc)}  has  the  linear  representation 

x(k)=  ^h(n)e(k-n) 

n~-~c«  5 

then 
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Bx{ C0!,C02)  = 

where 


oo 

//(©)=  XM«)^iVon 

n=-oo 

Example  2  -  Sum  of  independent  sinusoids.  If  4>  is  a  random  phase,  uniformly  distributed 
on  [0,27t],  then  the  bispectrum  of  the  process  (cos(XX-k|))}  is  identically  zero  for  any  radian 
frequency  X.  In  general,  if 

N- 1 

x(k)=  £  COS  (^„  +  <t>„)’ 

n  =  0 


where  { Xn }  is  a  given  set  of  frequencies  and  {(()„}  is  a  collection  of  i.i.d.  random  variables 
uniformly  distributed  on  [0,2k],  then  Bx( ©P©2)  =  0. 

Example  3  -  Phase  Coupling.  Consider  the  process 

x  (k)  =  cos  +  )+cos  ( kX2  +  $2)  +  cos  {kX2  +  (^3), 


where 


X2  <  Xj,  X3  =  Xj  +  X2,  <j>3  —  <|>|  +  <()2, 


and  4>i  and  <J>2  are  i.i.d.  random  variables  uniformly  distributed  on  [0,2k].  In  this  case,  the 
third-order  cumulant  sequence  is  given  by: 

C3(m,n)  =  -^{cos(X1n  +  X2m)  +  cos((X1  +  X2)n-X2m)  +  cos((X1  +X2)n  —  X]m^  J 
+ — jcos(X,m  +  X2n)  +  cos^X,  +  X2  )m  -  X2n)  +  cos^Xj  +X2)m-  X^jJ. 


Consequently,  the  bispectrum  of  the  process  evaluated  in  D 3  shows  an  impulse  at  the  point  (X^X^). 
This  is  in  sharp  contrast  to  the  previous  example  in  which  the  phases  were  not  coupled  and  the 
bispectrum  was  identically  zero.  It  is  for  this  reason  that  the  bispectrum  has  been  proposed  as  an 
efficient  means  of  detecting  quadratic  phase  coupling  in  a  process  [1]. 

Example  4  -  Bicoherence.  Let  {x(n) }  be  an  arbitrary  process,  and  consider  the  quantity: 

M0^’  “ 2 )  -  1  ,  —  —  x- 

This  quantity,  which  we  will  refer  to  as  the  bicoherence  of  the  process,  represents  a  normalized, 
scale  invariant  version  of  the  bispectrum  analogous  to  the  more  familiar  spectral  coherence 
function.  Unlike  the  spectral  coherence,  the  magnitude  of  the  bicoherence  of  a  process  is  not 
necessarily  bounded;  however,  it  can  be  shown  that  the  variance  of  a  conventional  estimate  of  the 
bicoherence  is  essentially  independent  of  the  power  of  the  process  .  Furthermore,  it  follows 
immediately  from  Property  3  above  that  the  magnitude  of  the  bicoherence  is  invariant  under  all 
linear  transformations,  and  the  phase  is  invariant  under  linear  transformations  with  linear  phase 


5 


response.  That  is,  if  {?(«)}  is  the  output  of  a  linear  system  with  frequency  response 

«(o>)  =  |//(co)|e'arg("<“)> 

driven  by  the  process  (x(n)},  then 

|py(cDi,  C02  )|  =  (a»i,  C02  )| , 

and  if  arg(#(co))  =  k  •  0),  where  k  is  an  arbitrary  constant  independent  of  (0,  then 

arg(py  ((Oj,  ©2 ))  =  arg(p *  (a)!,  to2)) . 

For  communication  purposes,  this  is  a  very  intriguing  property  of  the  bispectrum  (or  any  of  the 
higher-order  polyspectra)  since  it  implies  that  a  communication  scheme  based  on  amplitude 
modulation  of  the  bicoherence  will  be  insensitive  to  any  linear  distortions  introduced  by  the 
communication  channel.  Similarly,  a  communication  scheme  based  on  phase  modulation  of  the 
bicoherence  will  be  insensitive  to  distortions  with  linear  phase  response. 

A  Class  of  Polyspectrum  Estimates 

One  of  the  most  popular  conventional  polyspectrum  estimation  techniques  was  introduced  by 
Brillinger  and  Rosenblatt  in  [2],  To  define  the  Brillinger-Rosenblatt  (B-R)  estimate,  let  {xK(n)}  be 
a  data  sample  of  length  N  with  discrete  Fourier  transform  {XN(m) },  and  let  W  be  a  continuous, 
(&-l)-dimensional  frequency-domain  averaging  function  that  satisfies 

(i)  W(m)  >  0, 

(ii)  W(-c o)  =  W(w), 

(iii)  J  W(a)  da  =  1. 

R~l 

For  any  0  <  c  <  1,  define  the  scaled  version  of  W  by  Fl^^(co)  =  c^k  cp)j  Define  the 

generalized  periodogram  by: 

^N[nl’"'>nk-l)  =  ~j^^N{nl)^N{n2)'"^N{nk-l)^N{~nl  ~n2  nk-l)  , 

and  let  the  indicator  function  %  be  defined  by: 

0  if  =  0  (mod  2 k)  for  any  J  c:  1}, 

%(co)  = «  JeJ 

1  otherwise. 

Using  the  above  notation,  the  B-R  estimate  of  the  true  polyspectrum  S£  (co)  at  an  arbitrary  point  is 
given  by: 


6 


In  particular,  for  this  study,  we  always  used  an  averaging  parameter  value  of  c  =  0.5  and 
separable,  two-dimensional  rectangular  averaging  windows  of  the  form: 

w(q>1,<d2)=W1(o>1)w1(cd2), 

1  ..ii 
-  if  |C0|  <  7T, 

Wi((d)=|2jc 

0  otherwise. 

Hence,  all  of  our  bispectrum  estimates  took  the  particular  form: 


^w((0i,©2)  =  t|t  E  E 

Jy  n=  1  m=  1 


N-  1  N- 1 


2 7T77  1/9 

— yjt-N1'2^  Wx 
N112  1 


27tm  , ,, 
—T7J  — AT 1/2  (02 

A1/2  2 


Similarly,  all  of  our  bicoherence  estimates  took  the  form: 


0AT  (®1  ’^l)- 


2JV(m1,to2)  ’ 


where 


,  N-lN-l 

2»(®„0>2)  =  TT57tE  L  IV 

/V  rt  =  l  m=l 


2tt  /7  i/9 

- AT  2rr>  IV 

yy  1  12  W1  1 


27 T/72  i 

— nr- ^  ©2 

//* 2  2 


•|X„(n)XAf(m)X„(n  +  m)|. 

The  B-R  estimation  technique  was  employed  primarily  because  the  statistical  behavior  of  such 
estimates  is  well  understood.  Furthermore,  the  B-R  estimates  can  be  computed  very  efficiendy 
with  a  fast  algorithm  designed  and  implemented  by  the  principal  investigator  [3]. 


3 .  Signal  Design 

There  are  many  possible  signal  designs  that  one  might  use  in  order  to  encode  information  in  the 
bispectrum  or  bicoherence  of  the  transmitted  waveform.  In  essence,  each  possible  approach  can  be 
regarded  as  a  generalization  of  a  more  conventional  modulation  technique,  in  which  the  higher- 
order  moment  (HOM)  content  of  the  signal  is  available  in  addition  to  the  information  transmitted  in 
the  conventional  modulation.  In  situations  where  the  channel  performs  approximately  as  an 
additive  white  Gaussian  noise  (AWGN)  channel,  the  conventional  detection  scheme  is  often  nearly 
optimal  and  will  outperform  detectors  that  attempt  to  exploit  the  HOM  content  of  the  signal.  On  the 
other  hand,  in  situations  where  the  principal  source  of  channel  degradation  is  something  other  than 
AWGN,  such  as  fading  or  narrowband  interference,  considerable  performance  improvement  can 
be  realized  by  employing  a  detector  that  is  designed  to  exploit  the  higher-order  structure  of  the 
signal. 

In  this  study,  we  considered  baseband  signals  of  the  following  form: 


Where  {x, L  is  a  set  of  frequency  elements,  {<j);  }.=  {  a  set  of  phase  elements,  b  the  transmitted  bit 

value  ( 0  or  1),  and  mit)  a  spreading  signal,  which  will  be  discussed  in  detail  below.  We  will  also 
assume  the  following  relationships  among  the  various  elements: 

<]>!_(_ iy>  &  <|>2-(_iy>  uniformly  distributed  on  [2,rc], 

^'3-(-i)b  =  ^'i-(-i)1’"^^'2-(- iy” 

0  3  -  (-  1  f  ~  $  1  -  (- 1  f  ^  2  -  (- 1 1  ■ 

Notice  that  such  a  signal  is  really  a  generalization  of  both  a  frequency-shift-keying  (FSK) 
waveform,  in  which  three  phase-coupled  frequency  components  are  used  to  transmit  each  bit 
instead  of  a  single  component,  and  a  direct-sequence  spread-spectrum  (DSSS)  waveform,  in  which 

the  spreading  sequence  is  modulated  onto  a  multiple-component  carrier.  Similarly,  if  the  set  {A,t-} 
is  allowed  to  vary  randomly,  the  signal  becomes  a  generalization  of  a  frequency-hopped  spread- 
spectrum  (FHSS)  waveform.  For  expository  purposes,  we  will  generally  assume  that  the 
frequencies  are  fixed,  but  it  should  be  remembered  that  this  need  not  necessarily  be  the  case. 

We  will  be  concerned  primarily  with  the  sampled  version  of  this  signal,  denoted  by  the  sequence 
{*w(«)}„=0,  where 


xN  (n)  =  x  ( nT ) 

3 

=  m(nT)  ■  l  cos  (27tA,i_(_1)i,  nT+(. ) 

i=  1 
3 

&mN(n)  •£  cos(?if_(_iy, 

i=  1 

Here,  T  is  the  sampling  interval,  and  NT  is  the  length  of  the  bit  interval.  Notice  that  we  have 
implicitly  defined  the  spreading  sequence  {mN  (n)  }n=0 ,  and  the  set  of  normalized  frequencies 

if  ,  using  the  relationships: 

mN  (n)  =  m  (nT) , 

XN=2nTk. 

Since  we  will  almost  always  be  referring  to  the  digital  sequence  rather  than  the  analog  signal,  we 
will  usually  abuse  notation  somewhat  and  neglect  to  use  the  superscript  “AT  characters  except 
when  necessary  to  avoid  confusion.  Also,  we  will  generally  assume  that  a  bit  value  of  b  =  0  is 
transmitted. 

The  main  purpose  of  the  spreading  sequence  {m  (n ) }  is  to  increase  the  bispectral  bandwidth  of  the 
signal  in  order  to  improve  the  efficiency  of  the  bispectrum  estimates  employed  in  the  detector.  As 
discussed  above,  the  phase-coupled  carrier  itself  has  a  bispectrum  that  takes  the  form  of  an  impulse 
at  a  single  point  in  the  principal  domain.  This  turns  out  to  be  problematic  from  a  statistical 
standpoint,  at  least  when  using  B-R  estimates  of  the  bispectrum  and  bicoherence,  as  we  have  done. 
This  problem  can  be  alleviated  by  modulating  the  carrier  with  a  pseudo-random  sequence  having 
nonzero  bispectrum  and  bandwidth  chosen  to  match  the  “bandwidth”  of  the  bispectrum  estimator. 
In  particular,  if  the  bispectrum  estimates  are  computed  using  an  averaging  parameter  with  value 
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0  <  c  <  1,  then  the  bandwidth  of  the  spreading  sequence  in  the  digital  frequency  domain  should  be 
of  order  radians,  where  N  is  again  the  number  of  samples  in  a  bit  interval.  This  can  be 
accomplished  in  the  following  manner. 

Let  [h(k)}  be  the  impulse  response  of  a  real-valued,  low-pass  linear  filter  with  frequency  response 

H( to)  and  £  h2  (k )  =  1 ,  and  let  {hN(k)}  be  the  scaled  version  of  the  impulse  response  with 
*=-~ 

associated  frequency  response: 

1  —  c 

Hn((q)=N  2  //(iV1_cco),  -%Z(0<K. 

Let  {«(&)}  be  a  zero-mean,  i.i.d.  random  process  with  variance  a2,  third-order  moment  m3  *  0, 

and  bicoherence  (3  =  m3/cr^.  Filter  {u(k)}  with  {h/^k)}  to  produce  the  spreading  sequence 
{m(n)},  which  will  have  spectrum,  bispectrum,  and  bicoherence  that  satisfy: 

^(®)=a2K(co)f, 

■^7  if 

0  otherwise. 

If  {m{n)}  is  generated  in  this  fashion,  then  the  transmitted  signal  (x(n)}  has  spectrum,  bispectrum, 
and  bicoherence  that  satisfy  (assuming  b  =  0): 

Sx(Gi)=-^L^  +|tf;v(co+^)|  }> 

Bx((dl,d)2)^Y-HN((Ol-X1)HN((i)2-X2)HN(-(ii1-(ii2  +  X3),  near  (kltk2), 

|  Pi  (®1  >®2  )|  ="~3^  ’  near  (^1*^2)- 

<*s 

Now,  in  general,  the  bicoherence  of  the  signal,  when  considered  as  a  stationary  random  process, 
can  be  made  arbitrarily  large  by  choosing  an  i.i.d.  sequence  { u(k) }  with  arbitrarily  large  skewness 

(i.e.,  the  quantity  m3/a2)  as  the  input  to  the  filter  {/*/#)}.  However,  in  practice,  the  bispectrum 

of  the  signal  must  be  estimated  from  only  the  N  samples  available  during  a  single  bit  interval,  and  it 
is  not  difficult  to  show  that  the  sample  skewness  coefficient  of  a  data  sample  of  length  N  cannot  be 
larger  than  .  In  fact,  the  only  sequences  that  actually  attain  the  maximum  sample  skewness  are 
of  the  form: 


u  (k  )  —  {k )  yy  j, 


(1) 


where  {§N{k)  '  is  any  circular  shift  of  a  unit  impulse  sequence  of  length  N.  Since  we  have  the 

luxury  of  using  pseudo-random  rather  than  truly  random  sequences  in  our  modulation  scheme,  we 
will  always  choose  the  sequences  {u(k)}  to  be  of  form  (1). 

Finally,  for  the  purposes  of  this  study,  we  always  used  an  averaging  factor  of  c  =  0.5,  and  we 
employed  filters  {%(&)}  that  were  scaled  versions  of  a  canonical  low-pass  filter;  that  is, 

«(«)={' lf 

(0  otherwise. 

This  completely  describes  the  signals  used  in  this  study.  Henceforth,  we  will  refer  to  these  signals 
as  generalized,  phase-coupled  FSK  (GPCFSK)  waveforms  and  to  detectors  designed  to  exploit  the 
HOM  structure  of  such  signals  as  HOM  detectors. 


4 .  Performance  in  AWGN 

We  begin  our  analysis  of  the  performance  characteristics  of  bispectral  modulation  techniques  by 
considering  the  classical  AWGN  channel.  That  is,  we  assume  that  the  received  signal  sequence 
{ r(n) }  takes  the  form: 

r(n)=x  ( n)  +  w  (n), 

where  {x{n)}  is  a  GPCFSK  signal,  and  {w(n)}  is  a  zero-mean,  white  Gaussian  noise  process  with 
variance  o2w.  As  discussed  in  Section  2,  we  used  the  Brillinger-Rosenblatt  technique  to  generate 
estimates  of  bispectra  as  well  as  estimates  of  normalization  factors  that  were  in  turn  used  to 
compute  bicoherence  estimates.  Throughout  the  remainder  of  this  report,  we  will  make  frequent 

reference  to  the  asymptotic  distributions  of  the  estimates  |$ N  (X  l  ,X2 )|  ,  |$ N  (X4  ,X5  )|  , 
|2n(Xi,A,2)|2,  and  | 2n  (X4,A5)|2,  where  (hM)  is  the  normalized  frequency  pair  in  D3  associated 
with  the  bit  value  b  -  0,  and  (^4^5)  is  the  pair  associated  with  b  =  1.  If  we  let 

(*b,o=  Var  >^2)}’  m,o=  |e >^2) }|  5 

a5,l  =  Var  {^(^4^5)}’  1^5,1  =  |e  (^4  ’^5)}!  ’ 

^o=var{^(^i^2)}5  ^o=|e{^(A,i,X2)}|  > 

az,i  =  Var  {^(X,4,X,5)},  \izx  =  |e  \%n  (X,4,A,5)||  , 

then  it  follows  from  the  asymptotic  distribution  theory  developed  in  [2]  that  the  distributions  in 
question  can  be  approximated  as  follows: 


Here,  we  have  used  the  notation  (A, )  to  indicate  a  non-central  chi-squared  distribution  with  m 

degrees  of  freedom  and  non-centrality  parameter  X.  If  we  assume  that  separable  rectangular 
averaging  windows  and  an  averaging  parameter  of  c  =  0.5  are  used  in  all  B-R  estimates,  then  it 
follows  from  the  results  in  [2]  that, 
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where  represents  the  energy  transmitted  in  each  bit  sequence,  a3  represents  the  power  spectral 
density  (i.e.,  variance)  of  the  AWGN  process,  and  jt>  represents  the  signal-to-noise  ratio  (SNR), 
given  byyh  =  Eb/a2w. 


Now  consider  a  bicoherence  detector  that  employs  the  following  decision  rule  to  decide  whether 
b  =  0  or  b  =  1  has  been  transmitted: 


8({K«)0 


.  (Oif 

«)>)=  ^ 

1  if  0O2<03, 


where  we  have  adopted  the  notation: 

(X-i  ,A,2)|  ^2  jjjV_(^4>^S_)l 


K 

2»l 

kX.[,X2) 

Assuming  that  b  =  0  is  sent,  the  theoretical  probability  of  error  (Pe)  for  such  a  detector  is  given  by: 


where  the  symbol  Fmn  represents  an  F  distribution  with  m  and  n  degrees  of  freedom,  and  the  final 
expression  is  meant  to  be  interpreted  in  accordance  with  the  statement  “the  probability  that  a 
random  variable  with  distribution  F^i  is  less  than  a  constant  times  an  independent  random  variable 
with  distribution 

Similarly,  a  detector  based  on  the  bispectrum  rather  than  the  bicoherence  would  have  a  theoretical 
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Pe  given  by: 


«P 


To  get  some  idea  of  how  this  performance  compares  with  the  performance  of  more  conventional 
detectors,  we  note  that  if  one  ignores  the  phase  coupling  of  the  three  earners  comprising  the  signal, 
the  optimal  noncoherent  detection  strategy  is  to  compare  the  outputs  of  two  envelope  detectors 
matched  to  the  modulation  (i.e.,  spreading)  sequence  {m(k)}  and  appropriate  carrier  frequencies. 
The  performance  of  such  a  detector  (as  discussed,  for  example,  in  [4])  is  given  by: 


which  is  identical  to  the  performance  of  the  optimal  noncoherent  detector  for  conventional  binary 
FSK  signals  in  AWGN  [5].  The  theoretical  performance  of  the  bicoherence  and  bispectrum 
detectors  relative  to  this  optimal  detector  performance  is  illustrated  in  Figures  2  and  3  for  several 
values  of  the  parameter  N.  As  one  would  expect,  the  optimal  detector  outperforms  both  of  the 
HOM  detectors  regardless  of  the  length  of  the  bit  interval,  but  the  performance  of  the  bicoherence 
detectors  clearly  improves  as  N  increases. 
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Figure  2.  Theoretical  Performance  of  Bicoherence  Detector  on  AWGN 

Channel 
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Figure  2.  Theoretical  Performance  of  Bispectrum  Detector  on  AWGN 

Channel 


It  should  be  noted  that  the  bispectrum  detector  significantly  outperforms  the  bicoherence  detector  at 
higher  SNRs.  This  is  to  be  expected  on  AWGN  channels,  but  on  channels  that  exhibit  significant, 
frequency-selective  fading,  the  performance  of  the  bispectrum  detector  can  be  expected  to 
deteriorate  substantially  while  the  performance  of  the  Incoherence  detector  should  remain  relatively 
stable. 

In  order  to  validate  our  theoretical  predictions  of  the  performance  of  the  HOM  detectors  on  an 
AWGN  channel,  we  performed  several  simulation  studies.  The  results  of  the  simulations  are 
displayed  graphically,  relative  to  the  theoretical  performance  predictions,  in  Figures  4-11. 
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Figure  4.  Performance  of  Bicoherence  Detector  on  AWGN  Channel 

(N  =  256) 
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Figure  5.  Performance  of  Bicoherence  Detector  on  AWGN  Channel 

(AT  =  512) 
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Figure  6.  Performance  of  Bicoherence  Detector  on  AWGN  Channel 

(N  =  1024) 

(The  simulated  value  for  Pe  at  25  dB  is  an  approximate  upper  bound.  No  errors 
were  encountered  in  151,2u0  independent  trials.) 
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Figure  7.  Performance  of  Bicoherence  Detector  on  AWGN  Channel 

(N  =  2048) 
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Figure  8.  Performance  of  Bispectrum  Detector  on  AWGN  Channel 

(N  =  256 ) 
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Figure  9.  Performance  of  Bispectrum  Detector  on  AWGN  Channel 

(TV  =  512) 

(The  value  for  Pe  at  25  dB  is  an  approximate  upper  bound.  No  errors  were 
encountered  in  202,400  independent  trials.) 
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Figure  10.  Performance  of  Bispectrum  Detector  on  AWGN  Channels 

(N  =  1024) 

(The  simulated  value  for  Pe  at  25  dB  is  an  approximate  upper  bound.  No 
errors  were  encountered  in  15 1,200  independent  trials.) 
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Figure  11.  Performance  of  Bispectrum  Detector  on  AWGN  Channel 

(N  =  2048) 


As  these  figures  indicate,  the  correspondence  between  the  theoretical  performance  results  and  the 
simulation  results  is  similar  for  both  the  bicoherence  and  bispectrum  detectors.  In  particular,  for 
bit  intervals  of  length  1024  points  or  more,  the  theoretical  results  are  well  supported  by  the 
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simulations.  The  fact  that  larger  values  of  N  are  required  to  produce  simulation  performance  that 
agrees  with  theoretical  predictions  is  not  surprising  since  the  theoretical  results  are  all  based  on 
asymptotic  distribution  theory;  that  is,  the  theory  is  only  valid  as  N  goes  to  infinity. 


5 .  Performar.ee  in  Partial-Band  Jamming 

In  this  section,  we  consider  the  performance  of  bispectral  modulation  in  the  presence  of  partial- 
band  jamming.  We  restrict  attention  to  bicoherence  detection  strategies  since  they  are  somewhat 
easier  to  handle  analytically  than  bispectrum  detection  schemes,  both  from  a  theoretical  and 
simulation  standpoint.  In  any  case,  preliminary  analysis  indicates  that  the  performance 
characteristics  of  the  two  different  strategies  are  quite  similar. 

We  now  assume  that  the  received  signal  takes  the  form: 

r(n)=x{n)+j{n), 

where  [jin )}  is  a  zero-mean,  band-pass  interference  process  with  zero  bispectrum,  variance  G;  and 

bandwidth  equal  to  27ta  radians,  with  0  <  a  <  1.  This  implies  that  a  fraction  a  of  the  total 
available  bandwidth  is  being  jammed  during  each  bit  interval,  and  we  assume  that  the  jammer 
power  is  distributed  randomly  throughout  the  spectrum.  That  is,  in  any  bit  interval,  the  power 

spectral  density  of  the  jamming  process  is  identically  equal  to  G2/ a  on  its  region  of  support,  where 

the  region  of  support  varies  from  one  bit  interval  to  the  next  but  always  occupies  exactly  2na 
radians.  This  is,  of  course,  equivalent  to  the  assumption  that  the  jammer  remains  fixed  and  the 
user  randomly  varies  the  locations  of  the  signal  carrier  frequencies. 

Under  these  assumptions,  it  is  easy  to  see  that  the  total  probability  of  error  is  approximately: 
Pe=(l-a)3PeC  +  3a(l-a)2p,1  +  3a2(l-a)PeE  +  a3Pel3, 

where  Pei*  represents  the  probability  of  error  given  that  exactly  k  components  of  the  signal  have 
been  jammed  (i.e.,  a  fraction  kJ 3  of  the  signal  bandwidth  intersects  the  region  of  support  of  the 
jammer.)  Furthermore,  using  the  asymptotic  distribution  theory  developed  in  the  previous  section, 
it  can  be  shown  that: 
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Here  again,  Et,  is  the  energy  in  the  signal  during  one  bit  interval,  and  G2  is  the  total  jammer  power. 
The  worst-case  probability  of  error  for  the  bicoherence  detector  in  partial-band  jamming  occurs  at 
the  value  of  a  that  maximizes  the  expression  for  total  probability  of  error  given  above. 
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Unfortunately,  this  optimization  problem  does  not  admit  a  closed-form  solution,  but  an 
approximate  value  for  the  worst-case  probability  of  error  can  be  found  numerically  for  each 

combination  of  the  parameters  N  and  y^. 

As  a  comparison,  we  consider  the  performance  of  a  noncoherent  binary  FSK  (BFSK) 
communication  system  in  the  presence  of  partial-band  jamming.  It  can  be  shown  [5]  that  the 
probability  of  error  for  such  a  system  is  given  by: 


Similarly,  it  can  be  shown  that  the  probability  of  error  for  such  a  system  is  maximized  when 
a  =  a*,  where 


-  if  Y 
a 

l  if  yfr<4. 

Hence,  the  worst-case  probability  of  error  for  the  BFSK  system  is  given  by: 

—  if  Y**  4, 

yh 


K=< 


1  -y/4  . 


Te 


if  Y/,<  4. 


To  provide  a  simple  comparison  of  these  two  systems,  we  evaluated  numerically  the  approximate 
worst-case  probabilities  of  error  for  the  bicoherence  detector  for  several  different  combinations  of 

N  and  y t>.  The  results  of  this  evaluation  are  compared  graphically  with  the  worst-case  probabilities 
of  error  for  the  BFSK  system  in  Figure  12. 
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Figure  12.  Theoretical  Performance  of  Bicoherence  Detector  for  Worst-Case 

Partial-Band  Jamming 


As  the  figure  indicates,  the  theoretical  performance  of  the  bispectral  modulation  system  improves 
as  the  value  of  N  increases  but  stabilizes  around  N  =  1024.  Furthermore,  at  higher  values  of 
SNR,  the  bispectral  modulation  system  significantly  outperforms  the  BFSK  system  in  a  partial- 
band  jamming  environment. 

Once  again,  to  validate  the  theoretical  performance  predictions,  we  performed  several  simulation 
studies.  In  this  case,  we  simulated  the  performance  of  a  bicoherence  detector  operating  on  a 
GPCFSK  signal  in  the  presence  of  a  jammer  with  perfect  knowledge  of  the  signal.  That  is,  we 
assumed  that  the  bandwidth  of  the  jammer  was  exactly  matched  to  the  bandwidth  of  the  signal  for 
all  values  of  N  and  all  three  of  the  signal  components  were  certain  to  be  jammed.  Under  this 
scenario,  the  total  probability  of  error  is  exactly  equal  to  Pe|3,  which  we  use  as  the  baseline  for 
comparison  with  the  simulation  results.  The  results  of  the  simulations  are  presented  graphically, 
relative  the  theoretical  values  of  Pei3,  in  Figures  13-15. 
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Figure  13.  Performance  of  Bicoherence  Detector  Against  Omniscient  Jammer 

(N  =  256) 
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Figure  14.  Performance  of  Bicoherence  Detector  Against  Omniscient  Jammer 

(N  =  1024) 
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Figure  15.  Performance  of  Bicoherence  Detector  Against  Omniscient  Jammer 

(N  =  4096) 

Again,  as  the  figures  indicate,  for  bit  intervals  of  length  1024  samples  or  more,  the  theoretical 
results  are  well  supported  by  the  simulations.  Interestingly,  in  the  low  SNR  region,  where  the 
disparity  between  the  simulations  and  the  theoretical  predictions  is  the  greatest,  the  simulated 
performance  is  actually  considerably  better  than  predicted.  This  implies  that  even  for  low  SNR 
situations,  the  theoretical  results  provide  a  lower  bound  that  should  be  useful  for  performance 
predictions. 

6 .  Conclusions 

In  Phase  I  of  this  study,  we  have  completed  a  thorough  performance  analysis  of  one  possible 
design  for  a  digital  communication  system  based  on  bispectral  modulation.  We  have  analyzed  the 
performance  on  both  AWGN  channels  and  partial-band  jamming  channels  and  compared  that 
performance  to  the  known  performance  characteristics  of  more  conventional  modulation  schemes. 
As  expected,  the  performance  of  the  bispectral  modulation  technique  on  the  AWGN  channel  is 
inferior  to  the  performance  that  can  be  achieved  with  conventional  noncoherent  schemes,  but  the 
performance  in  the  presence  of  partial-band  jamming  is  significantly  better  than  conventional 
schemes  in  the  high  SNR  regimes.  This  implies  that  significantly  less  power  is  required  to 
maintain  a  high  data  rate  in  a  jamming  environment  when  using  bispectral  modulation  than  when 
using  conventional  modulation  such  as  BFSK.  Clearly,  this  has  implications  not  only  for  channels 
subject  to  intentional  jamming,  but  also  for  situations  where  two  or  more  frequency-hopped 
systems  are  operating  in  close  proximity  to  one  another,  or  where  two  adjacent  cells  in  a  cellular 
network  are  required  to  employ  some  frequency  reuse. 

Even  in  a  relatively  benign  AWGN  environment,  bispectral  modulation  may  offer  some  advantages 
over  conventional  approaches.  For  example,  if  covert  communication  is  desired,  a  GPCFSK 
system  could  be  employed  instead  of  a  simple  BFSK  system.  The  system  could  be  used  just  as  a 
conventional  BFSK  system  in  noncovert  mode;  that  is,  two  modulated  frequencies  would  be 
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assigned  to  each  of  several  users,  who  would  transmit  on  either  one  frequency  or  the  other 
depending  on  the  bit  value  they  wished  to  send.  In  covert  mode,  the  frequency  coupling  between  a 
subset  of  the  observable  users  would  be  coordinated  to  transmit  an  additional  data  stream  that 
would  only  be  detectable  with  an  HOM  detector.  The  enforced  coupling  among  the  chosen  subset 
of  users  would  cause  some  interference  with  the  remainder  of  the  users  but  such  interference 
would  normally  be  detectable  by  the  other  users,  who  would  simply  retransmit.  The 
nonparticipating  users  would  not  interfere  substantially  with  the  HOM  data  stream,  which  would 
essentially  be  undetectable  by  conventional  techniques. 

Similarly,  one  could  treat  a  GPCFSK  system  as  a  conventional  DSSS  system  in  which  each  user 
modulates  the  spreading  sequence  onto  multiple  carriers.  The  optimal  receiver  for  each  user  would 
then  incorporate  envelope  detectors  matched  to  the  spreading  sequence  at  each  of  the  appropriate 
frequencies.  Naturally,  such  a  detection  scheme  requires  that  the  detectors  acquire  and  track  the 
phase  of  the  chip  sequences  in  order  to  give  optimal  performance,  and  performance  could  be 
substantially  degraded  if  the  detectors  were  not  properly  synchronized.  In  contrast,  an  HOM 
detector,  which  could  be  run  in  parallel  to  the  conventional  detector,  requires  only  relatively  crude 
synchronization  in  order  to  maintain  stable  performance  characteristics.  In  this  case,  an  intelligent 
combination  of  the  two  detection  strategies  could  enhance  the  overall  performance  of  the  system. 

It  is  clear  from  the  results  of  this  study  that  polyspectral  modulation  techniques  in  general,  and 
bispectral  modulation  in  particular,  can  provide  improved  performance  over  conventional 
modulation  strategies  in  some  situations.  Considerably  more  research  involving  the  various 
modulation  techniques  available  will  be  required  to  determine  the  true  extent  of  the  performance 
advantages  and  the  feasibility  of  the  various  alternatives.  We  believe  that  our  Phase  I  results 
strongly  suggest  that  such  additional  research  is  warranted. 
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